arXiv: 1506.05191vl [cond-mat.str-el] 17Jun2015 


Origin of the glass-like dynamics in molecular metals k-(BEDT-TTF) 2X: implications 
from fluctuation spectroscopy and ab initio calculations 

Jens MiillerQ Benedikt Hartmann, Robert Rommel, Jens Brandenburg, and Stephen M. Winter 

Institute of Physics, Goethe-University Frankfurt, 60438 Frankfurt (M), SFB/TRfO, Germany 

John A. SchlueteJI] 

Argonne National Laboratory, Materials Science Division, Argonne, IL 60439, USA 

(Dated: June 18, 2015) 

We have studied the low-frequency dynamics of the charge carriers in different organic charge- 
transfer salts /^-(BEDT-TTF) 2 X with polymeric anions X by using resistance noise spectroscopy. 

Our aim is to investigate the structural, glass-like transition caused by the conformational degrees 
of freedom of the BEDT-TTF molecules’ terminal ethylene groups. Although of fundamental im¬ 
portance for studies of the electronic ground-state properties, the phenomenology of the glassy 
dynamics is only scarcely investigated and its origin is not understood. Our systematic studies of 
fluctuation spectroscopy of various different compounds reveal a universal, pronounced maximum 
in the resistance noise power spectral density related to the glass transition. The energy scale of 
this precess can be identified with the activation energy of the glass-like ethylene endgroup struc¬ 
tural dynamics as determined from thermodynamic and NMR measurements. For the first time for 
this class of ’plastic crystals’, we report a typical glassy property of the relaxation time, namely a 
Vogel-Fulcher-Tammann law, and are able to determine the degree of fragility of the glassy system. 

Supporting ab initio calculations provide an explanation for the origin and phenomenology of the 
glassy dynamics in different systems in terms of a simple two-level model, where the relevant energy 
scales are determined by the coupling of the ethylene endgroups to the anions. 


I. INTRODUCTION 

In the last years, the quasi-two-dimensional (quasi- 
2 D) organic charge-transfer salts (BEDT-TTF)2X, where 
BEDT-TTF (bis-ethylenedithio-tetrathiafnlvalene, com¬ 
monly abbreviated as FT) represents C6S8[C2H4]2, have 
attracted considerable attention as model systems for 
studying the physics of correlated electrons in reduced 
dimensions^®. The observed rich phenomenology of 
ground states is due to both the highly tunable na¬ 
ture of the correlation strength of the charge carriers 
and the strong coupling of the latter to intra- and in- 
termolecular vibrational modes of the underlying crys¬ 
tal lattice. The materials are highly anisotropic systems 
with alternating conducting and insulating layers, where 
the former contains the electron-donor FT molecules, 
for which the overlap of the molecular 7 r-orbitals of ad¬ 
jacent molecules leads to the formation of a quasi- 2 D 
conduction band. The /^-phase organic conductors k- 
(ET)2X with polymeric anions X“ are the most inten¬ 
sively studied systems of this family because of their su¬ 
perconducting, magnetic and dielectric properties, and 
in view of the current debate related to the Mott metal- 
insulator transition (MIT)P®. The latter can be accessed 
by varying the ratio of bandwidth W to effective onsite 
Coulomb repulsion U by means of physical or chemically- 
induced pressur^. For example, the ground state of the 
antiferromagnetic Mott insulator /^-(ET)2Cu[N(CN)2]Cl 
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(T/v = 27 K, denoted as n-Cl hereafter) can be turned 
into a superconducting state by either applying a mod¬ 
erate pressure of ^ 300 bar or by modifying the anion 
X: /i:-(ET)2Cu[N(CN)2]Br (k-Bt) and /^-(ET)2Cu(NCS)2 
(/^-CuNCS) are ambient-pressure superconductors with 
Tc = 11.6 K and 10.4 K, respectively. Likewise, stepwise 
deuteration, z.e., replacing hydrogen by deuterium in the 
terminal ethylene moieties [C2H4] of the FT molecules, 
drives the system tv-Bi from the superconducting towards 
the Mott insulating phas^^. In the past years, the in¬ 
trinsic structural disorder associated with the thermal 
dynamics of the terminal ethylene groups has been found 
to play an important role for the ground-state properties 
of the strongly correlated electrons. As shown schemati¬ 
cally in Fig.[^ the [C2H4] ethylene endgroups (EEC here¬ 
after) can adopt two different conformations, i.e., when 
viewed along the central C=C bond of the FT molecule, 
their orientation is either eclipsed (E) or staggered (S). 
For the /^-(ET)2X salts discussed in this work, the pop¬ 
ulation of the E and S states is thermally disordered 
at room temperature due to strong thermal vibrations. 
Upon cooling to low temperatures, the EEG tend to 
adopt one of the two possible conformations, depending 
on the anion and the crystal structure, ie., E for and 
K-C\ and k-Bt and S for /i:-CuNC^. However, in these 
three salts, the ordering of the EEG at low temperatures 
cannot be completed for kinetic reasons, since upon cool¬ 
ing through a characteristic temperature Tg their molec¬ 
ular motion slows down so rapidly that thermodynamic 
equilibrium cannot be reached. At this so-called glass¬ 
like transition, a short-range structural order becomes 
frozen-in. The glass temperature Tg depends in a char¬ 
acteristic way on the cooling rate qc, i.e., faster cooling 
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FIG. 1. (color online) (a), (b) ID coordination polymer an¬ 
ions Cu[N(CN) 2 ]Br and Cu(SCN) 2 , respectively, (c) Crystal 
structure of /^-(ET) 2 Cu[N(CN) 2 ]Br, side view showing the se¬ 
quence of conducting ET layers and insulating anion sheets, 
(d), (e) Schematic view of the relative orientations of the ET 
molecules’ ethylene endgroups (EEC) in the staggered (S) and 
eclipsed conformation (E), respectively. 


rates result in a hig her Tg and therefore a larger degree 
of quenched disordei^im^ Below this temperature, the 
thermal vibrations are frozen and no longer contribute to 
the specific heat and thermal expansiorP^. Thus, for the 
present materials, the glassy characteristics of the order¬ 
ing transition becomes evident in thermodynamic quan¬ 
tities, as frequency-dependent specific heafP^ and most 
pronounced in large anomalies in the a nisotropic, lin¬ 
ear coefficients of thermal expansiorPm^ exhibiting ther¬ 
mal hysteresis and under- and overshoot behavior in the 
warming curves. For thermal expansion measurements, 
a possible way to determine the glass-transition temper¬ 
ature Tg is to consider the intersection of warming and 
cooling curves taken at the same rates |g|, which results 
in —\q\ • dr/dT\Tg = 1, where r is the relaxation time of 
the corresponding molecular entities, see RefsP^^. This 
analysis yields = 72 K for k-C\ and = 75 K for k-Bt. 
k-CuNCS exhibits a more complex behavior, namely a se¬ 
quence of two transitions at Tg^ = 70K and Tg^ = 53 K. 
In contrast, no evidence of a distinct glass transition 
has been observed in other salts, such as the quantum 
spin liquid candidate /^-(ET) 2 Cu 2 (CN )3 (denoted as n- 
Cu 2 (CN) 3 )P^, despite thermal disorder of the EEGs in 
this material at high temperatures. In this article we 
aim to provide an explanation for these differences. 

From a fundamental point of view, the nature of the 
glass (or glass- lik^ transition is one of the most in¬ 
triguing and still unresolved problems in condensed mat¬ 
ter physic^i^M^. The phenomenon of glass formation 
occurs in a variety of different material^^^^^, whereas 
the tendency of building a glass former depends on the 
complexity of the contributing molecules. Besides ac 
specific heat, thermal expansivity and compressibility 
measurementJ^^M^, it is important to extract informa¬ 
tion of the frequency- and temperature-dependent dy¬ 
namics of the glass building processes. To that end, 
a number of sophisticated experimental techniques like 
viscosity measurements for liquidP^^ and broadband di¬ 


electric spectroscopy of solidP^ have been established. 
In particular the latter technique has become a pow¬ 
erful tool for studying the glassy dynamics of ’plastic 
crystals’ (or orientational glasses/^, a category to which 
the present materials belong, since they are character¬ 
ized by disorder with respect to the orientational degrees 
of freedom of the translationally ordered moleculeJ^. 
Dielectric spec trosco py has been successfully applied to 
many q uasi- and quasi-2D organic charge trans¬ 
fer salt d^ l ^^^ ^^, where both short-range relaxor-type as 
well as long-range ferroelectric order has been observed, 
mostly in charge-ordered and insulating phases. The 
glassy dynamics of the EEG rotational degrees of free¬ 
dom, however, has not been investigated by dielectric 
spectroscopy, since even the conductivity of semiconduct¬ 
ing K-C\ is too high for temperatures T > T^. 

In previous studies we were able to show that resis¬ 
tance 1//-noise spectroscopy provides a useful alterna¬ 
tive method for studying glassy molecular dynamics in 
conducting mate rials, where dielectric spectroscopy is 
difficult to appl^J^^^. A prerequisite is a sufficiently 
large coupling of the vibrational properties of the glass¬ 
forming entities to the electric conductivity and its fluc¬ 
tuations, which recently has been directly observed in 
NMR experiments^. The specific effects of this large 
coupling is twofold: 

(i) The anisotropic change of in-plane lattice parame¬ 
ters at Tjm results in a smaller average ratio tinter/^intra 
of the inter- and intra-dimer transfer integrals for more 
rapid cooling, which in turn leads to a smaller effective 
W/lM, i.e., the strength of electronic correlations can 
be controlled by varying the degree of frozen-in disor¬ 
der at Tg. In a recent work, we could demonstrate that, 
when thermally coupled to a low-temperature heat bath, 
a pulsed heating current through the sample causes a 
very fast relaxation with cooling rates at Tg of the or¬ 
der of several lOOOK/miiP^. The sudden freezing of the 
structural degrees of freedom causes a decrease of the 
electronic bandwidth W with increasing cooling rate, and 
hence a Mott metal-insulator transition for metallic sys¬ 
tems crossing the critical ratio (W/U) c of b andwidth to 
on-site Coulomb repulsion t/, see RefsP^l^. Due to the 
glassy character of the transition, the effect is persistent 
below Tg and can be reversibly repeated by melting the 
frozen configuration upon warming above Tg. Both by 
exploiting the characteristics of slowly-changing relax¬ 
ation times close to this temperature and by controlling 
the heating power, the m aterials can be fine-tuned across 
the Mott transit iorP^^^. 

(ii) The strong EEG vibrations lead to an enhanced 

sc atteri ng contribution at elevated temperatures above 
TJ 36 | 42 | their frozen-in configuration below Tg 

causes a random lattice potential, which gives rise to 
an additional contribution to the residual resistivity 
and affects the superconducting state and transition 
temperatur^^^®^. The influence of quenched disorder 
on the strongly-correlated 7r-carriers, which results in 
a Mott-Anderson transition, recently has become of in- 
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creasing theoreticaP^EH and experiment aP^H^ interest. 

Experimentally, disorder may be realized either by a 
systematic irradiation with X-ray^^ or by varying the 
cooling-rate qc = dT/dt through the glass-like transition 
at Tg. The latter type of disorder (in the EEG orienta¬ 
tions) is intrinsic in nature and can be controlled in a 
reversible way. Therefore, the cooling-rate dependence 
may be also utilized as a tool for studying the funda¬ 
mentals of the influence of randomness on the normal- 
and superconducting ground state of correlated-electron 
systems: cooling through Tg with increasing rate sys¬ 
tematically s uppresse s the superconducting Tc of k-Bt 
and In Ref. [121 we have described a 

simple model for tz-Cl and /^-Br, where the eclipsed and 
staggered EEG conformations are described as a double¬ 
well potential with an energy difference 2AE ^ 200 K 
between the E and S state and an activation energy 
barrier Ea of order 2400 — 3000 K. Using this model, 
we quantitatively estimate the amount of EEG units 
frozen in the non-equilibrium S orientation, as compared 
to the majority in the E state, for different cooling 
rates (see Eig. 3 in Ref. [ 12 ]). We find that for moder¬ 
ate and even the rapid cooling rates reported in the 
literature, the occupation of S conformation is below 5 %. 

Besides the well-established experimental phenomenol¬ 
ogy of the freezing-in of the [G 2 H 4 ] moieties and despite 
its significance for the electronic ground state properties, 
the glass-like transition has not been investigated in much 
detail and its origin in some /t:-(ET) 2 X salts with poly¬ 
meric anions - while no transition has been reported for 
other (ET) 2 X compounds - still remains puzzlin g. A lso, 
besides a recent NMR investigation on /^-CuNCS^^, the 
low-frequency dynamics of the molecular EEG motion 
itself has not been systematically studied. The fluctua¬ 
tion spectroscopy data presented in this work, due to the 
lack of dielectric spectroscopy at elevated temperatures, 
reveal for the first time characteristic properties specific 
for orientational glasses like a Vogel-Eulcher-Tammann 
behavior, which gives insight in the nature of the glass 
transition in organic charge-transfer salts and allows for 
a comparison with other glass-forming plastic crystals. 
Ab initio studies of the energetics of EEG rotation al¬ 
low the identification of several key distinctions between 
those materials displaying glassy freezing. 


iment are described elsewher^^. Care has been taken 
that spurious noise sources, in particular contact noise, 
do not contribute to the results. Resistance R and (re¬ 
sistance) noise power spectral density (PSD) Sr have 
been measured perpendicular to the conducting layers. 
Sr = Sv/P is calculated from the measured voltage 
noise PSD defined by 


5v(/) = 2 lim - 

T^oo 1 


T/2 

/ 

-T/2 




( 1 ) 


where I and SV{f) represent the current through and 
the fluctuating voltage drop across the sample, resepec- 
tively. Eor all samples, we have observed excess noise 
of general l/J'^-type, characterizing the intrinsic resis¬ 
tance (conductance) fluctuations. Typical spectra (raw 
data) which consist of 50 averages and a discussion of 
reproducibility and error bars are presented in a previ¬ 
ous publicatiorP^. The investigated samples have been 
continuously cooled down to 4.2 K with a slow cooling 
rate of about 3K/h. Eor /^-Dg-Br, we also applied a fast 
cooling rate of 300 K/h. Noise spectra have been taken 
while warming up the sample in discrete steps. The tem¬ 
perature stability during a noise measurement is better 
than 1 mK. The resistance values taken simultaneously 
with the noise measurements in discrete steps are found 
to be in very good agreement with the continuous mea¬ 
surements, see Ref. m 

Ab-initio calculations presented in section [TV| were per¬ 
formed using ORGAI^ at the B3LYP/def2-SV(P) level. 
Starting geometries for a series of ET salts were taken 
from room temperature crystal structures reported in 
Ref. [62l Galculations were performed on isolated ET 
dimers, with the interaction between the dimer and an¬ 
ion layer accounted for with a OPLS-aa forcefleld^ im¬ 
plemented in the qmmmopt option. Lennard-Jones pa¬ 
rameters were taken from the GROMAC^^ set, while 
the charge distribution of the polymeric anion layer was 
estimated from Mulliken analysis of B3LYP/def2-SV(P) 
calculations on small fragments. Eor each conformation, 
the coordinates of the EEG of interest were relaxed, along 
with the two bonded sulfur atoms, which typically dis¬ 
play large thermal ellipsoids at room temperature indica¬ 
tive of small disorder. Eor transition state calculations, 
ORCA’s hessian mode following algorithm was employed. 


II. METHODS 

Single crystals of plate- and rod-shaped morphology of 
/^-(ET) 2 X with X = Gu[N(CN) 2 ]C 1 , Cu[N(GN) 2 ]Br and 
Gu(NGS) 2 , as well as the fully deuterated variant /^-(Dg- 
ET) 2 Cu[N(CN) 2 ]Br (denoted as /^-Dg-Br) were grown 
by electrochemical crystallization according to literature 
methodPH^. 

Low-frequency fluctuation spectroscopy measurements 
have been performed in a five-terminal setup using a stan¬ 
dard bridge-circuit ac techniqu^^. Details of the exper¬ 


III. EXPERIMENTAL RESULTS 

III.l. Signature of glass-like structural dynamics in 
resistance noise spectroscopy — the DDH model 

Eigurej^ shows the noise PSD of the resistance fluctua¬ 
tions, Sr/R^^ taken at 1 Hz and normalized to the value 
at T = 250 K, for five compounds of the /^-(ET) 2 X family 
located at different positions in the generalized phase dia¬ 
gram, and thus having different electronic ground states. 








4 



FIG. 2. (color online) Resistance noise PSD, Sr/R^, taken 
at 1 Hz and normalized to the value at 250 K, vs. T for five k- 
(ET) 2 X salts located at different positions in the generalized 
phase diagram: between insulating r-C\ and metallic /^-Hs-Br 
and ac-CuNCS, pressurized /^-CP and fully-deuterated K-Dg- 
Br are located in the critical region close to the Mott MIT. For 
A^-Hs-Br data for two different samples are shown. Curves are 
shifted for clarity. Data have been taken in discrete steps in 
warm-up measurements after continuously cooling down the 
samples with slow cooling rates. 


K-C\ is insulating under ambient conditions and is shifted 
to the critical region close to the Mott MIT by appli¬ 
cation of physical pressure (/^-Cl*, same sample as in 
Ref. [36]). Likewise, chemically-induced pressure in fully- 
deuterated hz-Dg-Bi (same sample as in Ref. |42|) drives 
the ambient-pressure superconductor /^-Hg-Br close to 
the Mott transition. /^-Hg-Br and a^-CuNCS are metallic. 
Strikingly, in spite of the very different electronic ground- 
state properties determined by the electronic correlation 
strength WjU^ the behavior at elevated temperatures of 
the different compounds is very similar, namely governed 
by a pronounced maximum in the noise power centered 
around T ^ 100 K. Besides the broad maximum, the 
data sets display an anomalous feature at the thermo¬ 
dynamic glass-transition temperature T^, i.e. a peak or 
small shoulder in Sr/R^{T). It is natural to conclude 
that the overall temperature dependence is due to the 
structural similarities of the different compounds. 

Besides the magnitude of the 1 //-type noise shown in 
Fig. ID it is interesting to analyze its frequency depen¬ 
dence. In a previous stud}!^ we have demonstrated that 
the temperature dependence of the observed excess noise 
can be explained by the phenomenological Dutta-Dimon- 
Horn (DDH) modeP^, which assumes that the l//°^-type 


noise is caused by many independent ‘fluctuators’ - each 
contributing a Lorentzian PSD - acting in concert. These 
fluctuators may be identified with thermally-activated 
two-level processes characterized by a relaxation time 
r = To exp (F^//cbT), where tq ~ (10 “^^ — 10“^^) s is 
an attempt time in the order of the inverse phonon fre¬ 
quencies, E the activation energy of the process and /cp 
Boltzmann’s constant. The voltage autocorrelation func¬ 
tion for such a two-level system, 

Av{t) = {SV{t + T) -SVit)) = (W)2exp(-|t|/r), (2) 


is purely exponential, while the PSD, according to the 
Wiener- Khintc hine theorem, is a Lorentzian function of 
frequenc}!^^® 

_ 4r 

Svif) = 4 / Av(t) cos {2Trft)dt = {5VY ^ 4 ,^ 2 /^^ ' 

(3) 

In a particular ’noise window’ the action of an individual 
fluctuator may be enhanced over the 1// background, 
such that the kinetics and characteristic energy of this 
process can be determinecP^. 

The DDH model assumes that the a priori unspecified 
fluctuators linearly couple to the measured resistance 
fluctuations. The main ingredient of the model is a distri¬ 
bution of activation energies D{E)^ which governs both 
the strong temperature dependence of the magnitude of 
the noise, Sr/ R^{f^T)^ and any deviations from per¬ 
fect {a = 1) 1/f behavior, i.e. the temperature depen¬ 
dence of the l//^-noise frequency exponent, a{T). In 
this model, it is 



and 


Tq exp {E/k^T) 

1 + Tq exp (2E/kR,T)ATT‘^p 


D{E)dE, 

( 4 ) 


<^calc(/,^) = 1 


1 / ain|i(/,r) 

ln27r/ro v 91nT 


d\ng{T) 

ainT 


Assuming that D{E) is a slowly varying function 
pared to thermal energies it follows that 


D{E) cx 


2'Kf 1 Sr 
kBTg{T) i?2 


if,T). 



( 5 ) 


com- 


(6) 


Due to the large logarithmic factor in = 

—fepT In (27r/ro), ordinary activation energies in solids 
can be accessed^. The strength of this approach is that 
- if D{E) in Eq. @ does not explicitly depend on tem¬ 
perature - the assumptions of the DDH model can be 
checked for consistency by comparing the measured fre¬ 
quency exponent 


s{T) = - 


ain|f(/,T) 


dlnf 


( 7 ) 


with the predictions <acaic(^) of fho model from Eq. 
We discuss here a generalization of the DDH modeP 
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which accounts for possible deviations of the predicted 
(calculated) frequency exponents from the measured 
data: a function g{T) may account for an explicit temper¬ 
ature dependence of the distribution of activation ener¬ 
gies. The physical meaning is that the coupling constant 
between the random processes and t he resistance, and 
hence the total noise magnitude {dRY = 5'i^(/)d/, is 

not independent of temperature. In the simplest case of a 
constant offset (vertical shift) of (acaic(^) as compared to 
o^meas{T)^ g{T) satisfies a simple power law g(T) = aT^^ 
where the exponent b is given by the vertical shift and 
describes the temperature dependence of the number 
and/or coupling strength of the fluctuators. Fig. i(a) 
shows, exemplarily for /^-Dg-Br, an excellent agreement 
of the measured data with the DDH modeP^, when for 
higher temperatures T > Tg the almost constant devia¬ 
tion between the measured and calculated frequency ex¬ 
ponent is accounted for by a power law g{T) = aT^^ 
where b = —3/2 gives a reasonably good description. 
Hints for a physical interpretation of this power-law be- 
vaviour comes from the temperature dependence of the 
resistance shown in the inset, where above Tg a scattering 
contribution oc can be identified, suggesting that 

the responsible mechanism may be linked to the function 
g{T) via the number and/or strength of the fluctuators, 
similar as has been argued for the contribution of spin 
fluctuations in certain manganiteJ^. In our case, the 
fluctuators at elevated temperatures can be associated 
with structural excitations involving the EEG, the cou¬ 
pling of which to the resistivit}!^ varies as see 

Ref. [721 

The non-monotonic temperature dependence a{T) 
shown in Eig.[^(a) is related to the energetic fingerprints 
of the fluctuating entities, since for a greater and smaller 
than 1, dD{E)/dE > 0 and dD{E)/dE < 0, respec¬ 
tively. The good agreement of Q^meas(^) with the gen¬ 
eralized DDH model allows for extracting D{E) from 
the temperature dependence of the magnitude of the re¬ 
sistance noise PSD by using Eq. Eigurej^b) shows 
D(E, / = 1 Hz), calculated with g = for T >Tg. 

The features below E ^ 0.1 eV (see arrow) corresponding 
to the sharp peak in a{T), which is also accounted for 
by the DDH model, are due to a sudden slowing down 
of the charge carrier dynamics in the vicinity of th e crit- 
ical endpoint of the Mott MIT as discussed in Refs.S^^^. 
Besides these electronic correlation effects at lower ener¬ 
gies, D(E, / = 1 Hz can be reasonably well described by a 
small linear term and a peak function (a Lorentzian func¬ 
tion is shown in Eig.j^b); a Gaussian peak function works 
equally well around the maximum, see Eig.[^ below) with 
the maximum at Emax ~ 262 meV (~ 3040 K) and full- 
width-at-half-maximum of ^ 160 meV (^ 1860 As 

we have discussed in our previous studies, Emax corre¬ 
sponds to the activation energy Ea of the orientational 
degrees of freedom of the ET molecules’ terminal ethy¬ 
lene groups, as deter mined e.g. by NMR, specific heat 
and thermal expansio d^^ 1 14 | 74| 7^^ Thus, in this case the 
relevant fluctuators can be identified a posteriori as the 



FIG. 3. (color online) (a) Temperature dependence of the 
l//“-noise frequency exponent, o(T), for n-Dg-Bv. Blue dots 
represent the data, the colored curves show the prediction 
Ocaic after Eq. 0 with g{T) — const, (solid red line) and for 
the generalized DDH model with g{T) = aT^ for T > Tg 
with b = —3/2 (dashed red line), b = —1 (dark grey) and 
b = —2 (light grey). Inset shows the temperature dependence 
of the resistance, which above Tg has been fitted to R{T) = 
RoTRiT-\-R 3 / 2 T~^^‘^ . (b) Distribution of activation energies 
D(E,f = 1 Hz) in the generalized DDH model, i.e. using 
g{T) = aT~^^‘^ for T > Tg in Eq. Line is a fit to a 
Lorentz peak function and a linear term and serves as guide 
to the eye (a Gaussian fit works equally well in the region 
around the maximum, see also Fig.[^. Thin arrow points 
to enhanced fluctuations induced by electronic correlation^^. 
Thick arrow marks the Energie Emax ~ 260 meV. 


ET molecules’ EEGP^. Thus, the strong resistance fluc¬ 
tuations at around 100 K shown in Eig.[^ originate in the 
strong coupling of the electronic degrees of freedom to 
the glass-like structural excitations, see alscP^. The cor¬ 
responding change in the frequency exponent a{T) upon 
cooling from a < 1 to a > 1 corresponds to a strong shift 
of spectral weight to lower frequencies. This is related to 
the expected slowing down of molecular dynamics when 
approaching the glass transition and consistent with our 
interpretation discussed below, namely that the observed 
charge carrier dynamics is caused by the slow, so-called 
Qf-process of the glass-forming EEG rotations embed¬ 
ded in the anion structure forming a ’ cag e’ environment 
through short G—H • • • anion contacts^. In dielectric 
spectroscopy, the ^-relaxation peak of the frequency- 
dependent dielectric loss - which may be connected to 
^(/) by a fluctuation-dissipation relation - shifts with 
decreasing temperature to lower frequencies in a similar 
way. 


III.2. Analysis of the dynamic transport properties 

Next we show that noise measurements allow for study¬ 
ing characteristic properties of the glassy system as, e.^., 
a possible non-Arrhenius behavior, which so far has not 
been investigated for the present materials. The fre¬ 
quency range of the 1 //-type fluctuations for the present 
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materials is naturally limited to about 1 mHz — 1 kHz. 
In this range, the relaxation is dominated by the a- 
relaxation, which is the slowest dynamic process in glass 
formers. In contrast to other, faster processes of differ¬ 
ent origin (as, e.^., the high-frequency /3-process, which 
can be visualized as a ’rattling’ movement of a particle 
in the transient ’cage’ formed by its surrounding), the 
(^-relaxation is associated with the forming and decay of 
the cage, involving the cooperative movement of many 
particle^, in our case the ET molecules’ EEG and their 
anion environment. The increase of the (^-relaxation time 
during cooling and deviations from a thermally-activated 
behavior in glass-forming materials is often parameter¬ 
ized by a n emp irical law employed by Vogel, Eulcher and 
TammanrP^l^ 


T = To exp 


D • Tvft 
T — Tvft 


(8) 


where tq is a prefactor, D is the so-called strength pa¬ 
rameter and Tvft is the Vogel-Eulcher-Tammann tem¬ 
perature. The divergence of the relaxation time at Tvft 
cannot be observed in a real experiment because the sys¬ 
tem falls out of equilibrium when the relaxation time 
becomes larger than the time scale of the experiment. 
The model suggests a ’hidden’ phase transition below 
the glass transition at T^, which, however, is controver¬ 
sially discussed. The strength parameter D can be used 
to classify the glass forming materials in terms of the 
deviations of t(T) from a t hermally activated behavior, 
for which Tvft = So-called ’fragile’ glass for¬ 

mers reveal small values of D (typically D < 10) with 
large deviations from an Arrhenius behavior, whereas 
’strong’ glass formers {D > 10) follow an Arrhenius 
law more closely. This classification scheme has proven 
very useful since many properties of glass-forming ma¬ 
terials are correlated with their strength or fragility. 
Typical fragile glass formers are substances with nondi- 
rectional interatomic/intermolecular bonds, e.g. molten 
salts. Strong systems, i.e. those which show a strong 
’resistance’ against structural degradation when heated 
through their supercooled regime, are found among the 
non-hydrogen bonded network melt^^. Another conve¬ 
nient measure for the strength of the fragility can be ob¬ 
tained by focusing o n the l ong relaxation time end of the 
strong-fragile patterrP^EH. The (dimensionless) slope at 
Tg is given by 


m = 


dlogr 


d(Tg/T) 


T=Ta 


(9) 


and is related to the strength parameter D by 

m = (D/lnlO) • (TvFx/Tg) • (1 - TvFx/Tg)-' (10) 
= TOmin + TO^i„lnlO/D 

where rrimin ~ 16 is found to be the lower limit of 
fragilit}E21. Eor a large class of glass-forming liqui^ the 
upper bound for the fragility index is rumax = 17Cp^. 
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FIG. 4. (color online) Data (circles) of the normalized re¬ 
sistance noise PSD recalculated from the l//“ frequency ex¬ 
ponent o(T) shown in Fig.|^a) and Gaussian fits (lines) vs. 
temperature for various frequency of Ac-Dg-Br. (The distribu¬ 
tion of activation energies D(E) in the DDH model shown in 
Fig[^b) are derived from the data for / = 1 Hz.) The black 
squares connecting the maxima of the fits represent the tem¬ 
perature dependence of the relaxation time characteristic for 
the EEG structural excitations. 


Here, we employ the Vogel-Eulcher-Tammann expres¬ 
sion, Eq. (§, in order to describe the temperature depen¬ 
dence of the peak in the l//^-noise PSD caused by the 
glass-like structural excitations at a certain frequency. 
The frequency dependence is essentially given by the 
noise exponent (t(T) shown in Eig.l^a) for /^-Dg-Br. The 
recalculated noise amplitude Sr/^{T) for different fre¬ 
quencies is shown in Eig.[^ The black squares that con¬ 
nect the maxima of the Gaussian fits thus represent the 
temperature dependence of the relaxation time charac¬ 
teristic for the EEG • • • anion structural excitations. In 
Eig.j^a) we plot these maxima at different frequencies in 
an Arrhenius plot for different systems. Strikingly, the 
curves follow a thermally-activated behavior for high fre¬ 
quencies, in agreement with earlier measurements of ther¬ 
modynamic propertieJ^. Eor lower frequencies, however, 
a clear deviation from an Arrhenius behavior is observed 
and the curvature can be described well by the Vogel- 
Eulcher-Tammann law, Eq. The corresponding slow¬ 
ing down of the motion of the structural units being sig¬ 
nificantly stronger than predicted by the Arrhenius law 
is a typical characteristic of a glassy systerrP^. As we 
will discuss below, a possible explanation for this behav¬ 
ior is a growing number of correlated molecules moving 
cooperatively which leads to an increase of the apparent 
activation energy for decreasing temperatures. 

We note that the definition of the glass transition temper¬ 
ature is not unambiguous since it depends on the cooling 
rate. A possible procedure for defining Tg from ther¬ 
modynamic measurements is described above. In spec- 
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troscopy measurements, Tg is often defined as the tem¬ 
perature wherethe relaxation time becomes r = 100 s, 
see Equivalently, in order to have a well- 

defined transition temperature for different materials, Tg 
is determined as the temperature, where the system falls 
out of equilibrium for q = lOK/min. Although we have 
employed cooling rates different than q = lOK/min, the 
corresponding frequency / = lA27rr) ~ 1.6 mHz is 
shown as a horizontal line in Fig.l^a). As expected, the 
Vogel-Fulcher-Tammann temperature Tvft is lower than 
the glass transition temperature Tg with the VFT-fits in 
Fig. E a) for the various slowly-cooled systems yielding 
Tvft = 65K for k-CI, Tvft = 64K and 85 K for slowly 
{q = 3K/h) and fast {q = 300 K/h) cooled /^-Dg-Br, re¬ 
spectively, and 68 K and 64 K for two different samples 
of /^-Hg-Br. The cooling-rate dependence of the VFT be¬ 
havior, observed for tv-Dg-Bi, is unusual since for classical 
glass-forming systems the curves for different q should co¬ 
incide for T > Tg. Rather than being a purely dynamical 
effect, we suggest that for the present glass-like materi¬ 
als different cooling rates produce different glass-forming 
molecular environments (see discussion below) resulting 
in different values of Tg. 

Orientational glasses usually are ’strong’ glassformers, in 
the sense that they exhibit approximately an Arrhenius 
behavior in a plot of the viscosity rj of the liquid In 77 vs. 
TgjT^ indicative of a temperature-independent activation 
energ}!^. Fragile liquids on the other hand exhibit super- 
Arrhenius behavior, i.e. their effective activation energy 
increases as temperature decreases. From our analysis 
of the noise measurements, we find a wide range of the 
strength parameter with T = 5 for k-C\, D = 8.6 and 
D = l.S for slow- and fast-cooled /^-Dg-Br, respectively 
and D = 2.S and 5.2 for two different samples of /^-Hg-Br 
(only one shown in Fig.[^. Apart from obvious sample- 
to-sample dependencies and errors related to the limited 
number of data points, values below D = 10, consis¬ 
tently observed for a total of 6 samples, reveal that the 
present charge-transfer salts are rather fragile glasses, in¬ 
dicating significant cooperativity between EEG degrees 
of freedom. Accordingly, the values for the slope are de¬ 
termined to be in the range of m ^ 100 — 200 for the 
different samples, i.e. in each case much larger than the 
lower limit of fragility corresponding to m = 16 (Ref. [22]). 
For comparison, in Fig.j^b) we plot the data for /^-Hg-Br 
and K-C\ in the classification scheme with other plastic 
crystals or orientational glassed. 

Besides the VFT behavior, a typical characteristic of 
glass formers is non-exponential relaxation, which may be 
caused by a distribution of relaxation times and usually is 
seen in an asymmetric broadening of the peaks in the di¬ 
electric loss^. In general, the fluctuation-dissipation the¬ 
orem relates the noise PSD of a quantity x to the imag¬ 
inary part of the complex susceptibility where 

the real function A{t) is the response function of the 
systenP^ 

S,{f) = 2k^T^^. ( 11 ) 

UJ 



FIG. 5. (color online) (a) Arrhenius plot of the noise peak fre¬ 
quency as shown exemplarily in Fig. I^for selected systems of 
/^-(ET) 2 X (see legend). Light grey lines are fits of thermally- 
activated (Arrhenius) behavior to the high-frequency data, 
whereas dark grey lines are VFT fits after Eq. ^ to the low- 
frequency data. The horizontal line equals r = 100 s. (b) The 
data for /^-Hg-Br and k-C\ in comparison to various other 
orientational glasses (taken from Ref. I83|) . 


Despite our limited frequency range, transforming the 
data of Fig. [^according to Eq. ( pT| reveals - as expected 
for a glass-forming system - asymmetric and considerably 
broader peaks than 1.14 decades half width given by the 
Debye equatiorP^. 


IV. DISCUSSION 

After having described the phenomenology of the glass¬ 
like EEG ordering in resistance noise spectroscopy, we 
now discuss the question of the occurrence and origin of 
the effect in /i:-(ET) 2 X. While enhanced electron scat¬ 
tering associated with thermally-excited structural fluc¬ 
tuations in the EEG orientational degrees of freedom 
may be a common feature of many (ET) 2 X salts, a pro¬ 
nounced glass transition has been reported in the lit¬ 
erature for the present /^-Cl, /^-Br, and /^-GuNCS sys¬ 
tems. In all three materials (and deuterated variants), 
the planar anion layer consists of one-dimensional Gu(I) 
coordination polymers (repeat unit X = Gu+[Y“][Z“]) 
with bridging ligands (Y) completing the backbone, and 
short side-chains formed by terminal ligands (Z) coor¬ 
dinated to each Gu(I). For example, the isostructural 
K.-CI and hz-Bi crystallize in the space group Pnma^ 
with zigzag anion chains consisting of Y = dicyanamide, 
[N(GN) 2 ]“, and Z = Cl“ or Br“, while k-CuNCS (for 
which Y = Z = thiocyanate, (SGN)“) crystallizes in the 
P2i space group. Other salts in this latter space group, 
such as /^-(ET) 2 Cu(GN)[N(GN) 2 ], i.e. Y = (GN)-, Z = 
[N(GN) 2 ]“, apparently show no evidence for structural 
disorder in the EEGs, let alone a glass transition. Gu- 
riously, thermal disorder of the EEGs has been found 
at room temperature in the P2i/c quantum spin liq¬ 
uid material /^-Gu 2 CN 3 , although thermodynamic mea¬ 
surements suggest a smooth ordering into the staggered 








conformation at low temperatures without a pronounced 
glass transition!^ In this latter material, the anion layer 
forms a more rigid 2D coordination polymer in which 
(CN)“ anions bridge adjacent -CuCN- chains. Given 
the structural similarities of the ET layers in each ma¬ 
terial, the potential for glass formation must originate 
in the specific nature of the coupling between the anion 
layer and EEG degrees of freedom. Indeed it has been 
suggested that such coupling is necessary for the strong 
cooperative dynamics between EEGs indicated by the 
here observed fragile glass response (D < 10), as the an¬ 
ion layer is nec essary to mediate EEG-EEG interactions 
over long range^®^!^. In order to explore the quali¬ 
tative effects of the EEG-anion coupling between these 
materials, the energetics of the E and S conformations 
has been investigated at the B3LYP/def2-SV(P) level us¬ 
ing ORGA, and treating the anion-cation interaction in 
the static OPLS-aa approximation, as described in Sec. 
[m The results are summarized in Eig. [^c) in terms of 
the computed ground states (E or S), energy differences 
2 AE = IE’s — EeI, and activation energies estimated 
from {Ea = Ets ~ (^s + ^e)/ 2), where ‘TS’ refers to 
the transition state. In all cases, the computed ground 
state agrees qualitatively with the observed conforma¬ 
tional preferences. 

Eor the isostructural k-C\ and k-Bt, the EEGs exhibit¬ 
ing disorder are all crystallographically equivalent, and 
were computationally found to prefer the E conformation 
by 2AE/kB ^ 750 K and 520 K per ET molecule for k- 
G1 and /^-Br, respectively. Despite qualitative agreement 
with experiment, these calculated energy differences are 
somewhat larger than the value of 2 AE//cb ^ 210 K es¬ 
timated from the temperature dependence of the relative 
E/S occupancy of k-Bi for fixed cooling rat^^^. The com¬ 
puted Ea values of 1930 K (Gl) and 2200 K (Br) agree 
well with the experimental range of Ea/kB = 2000—3100 
K as determined e.g. by NMR, specific heat, thermal ex¬ 
pansion, and resistivity measureme nts of str uctural re¬ 
laxation kinetics on these material^m^^^TlIZll por both 
salts, more than half of the computed energy differences 
2AE result from strong van der Waals coupling between 
the EEG and anion layers, which is dominated by close Z 
• • • H (or D) contacts between each EEG and two terminal 
ligands in adjacent anionic chains (Eig. [^a)). In turn, 
each Z = Br“ or Gl“ interacts with four EEGs, two in 
each adjacent ET layer. These terminal ligands exhibit 
significant thermal motion at room temperature, partic¬ 
ularly parallel to the chain direction, as evidenced by a 
large root mean square displacement of ^/Uaa = 0-29 A 
for /^-Br and 0.25 A for k-C\. The coupling of this motion 
with the EEG rotation in nearby ET molecules provides 
a possible mechanism for the cooperative behavior im¬ 
plied by the observed deviations from the Arrhenius law 
in Eig.|^ In this way, it is the collective motion of EEG 
and anions that freezes out at Tg instead of being simply 
a property of the vibrational degrees of freedom of the 
EEG within each individual, separated ET layer. Re¬ 
cently, it has been shown that for glass-forming liquids 



Computed p 
Ground State: ^ 


S/S 


S E/S 


FIG. 6. (color online) (a) Optimized geometries of the termi¬ 
nal EEG for k-Bt in the staggered (S), eclipsed (E), and in¬ 
termediate transition state (TS) conformation, (b) Schematic 
two-level potential energy surface defining Ea and AE. (c) 
Computed ground state, Ea and 2AE values for selected 
(ET) 2 X salts. For those salts with nonequivalent EEGs, the 
specific sublattice is indicated by an ” A” or ” B”. Only those 
salts on the left of the figure, indicated by a green background, 
display glass transitions. 


an increase of the number of correlated molecules causes 
the noncanonical (non-Arrhenius) slowing down of the 
molecular motion upon approaching which we also 

observe for the present materials. Eurt her more, the ac¬ 
tivation energy was found to be approximately propor¬ 
tional to the number of cooperatively moving structural 
units. This may also be the case for the EEG motion 
and may be related to the change of the lattice constant 
with temperature. Spatially-resolved measurements, as 
e.g. scanning tunneling microscopy, would be highly de¬ 
sirable in order to investigate if there is an inhomoge¬ 
neous distribution of the EEG orientations with ordered 
domains as suggested in Ref. [45l 

The computational results for hz-Cl and hz-Bi may be 
contrasted with those of the ’spin liquid’ compound k- 
GU 2 GN 3 . In the latter material, all disordered EEGs are 
also crystallographically equivalent, due to the presence 
of inversion centres in the P2i/c space group. However, 
appearing on such an inversion centre in the anion layer 
is a (GN)“ ion whose orientation (-GN- or -NG-) is statis¬ 
tically disordered. Eor this reason, the EEG-anion cou¬ 
pling was computed from the average of both cyanide 
orientations. In agreement with experimental findings, 
a preference for the S conformation was found, albeit 
with a much larger 2AE/kB ^ 2410 K, suggesting re¬ 
duced metastability of the higher energy E conformation. 
Taken together with the absence of soft terminal ligands 
in the anion layer to mediate cooperative motion, this 
finding warrants against glassy freezing in this material. 

Eor /^-GuNGS and /^-(GN) 2 N(GN) 2 , the lower symme¬ 
try of the P 2 i space group results in two inequivalent 
EEGs with the potential for thermal disorder. As shown 
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FIG. 7. (color online) Preferred conformations of the EEGs 
for (a) k-Bt and (b) k-CuNCS in relation to the nearby anion 
layer. Dashed lines indicate close EEG • • • Z contacts, while 
blue arrows indicate soft vibrational degrees of freedom of the 
terminal Z ligands suggested to couple to the EEG rotation. 
In Av-Br, all EEG are crystallographically equivalent, while 
/^-GuNGS consists of two unique A- and B-sublattices. 


of freedom is suggested by significant structural changes 
in the anion layer in the vicinity of 100 K and 60 K in 
k-CuNCS, which preempt the glass transitions at Tg^ = 
70 K and Tg^ = 53 K, respectivel}i^. These changes are 
particularly dramatic in modulation of the average Cu-N 
distance presumed to be due to motion of the terminal Z 
= SCN group. 

It has been previously suggestecP^^^ that the arrange¬ 
ment of the ET molecules and conformational preferences 
in /i:-(ET) 2 X salts with polymeric anions result from a 
particular pattern of empty spaces in the anion layers, 
allowing the donor molecules to fit in between, thus min¬ 
imizing the van-der-Waals interactions. From this per¬ 
spective, glassy freezing of the EEG rotation requires 
that the EEG-anion interaction is relatively ineffective 
at distinguishing the S and E states, allowing for similar 
energies of both states, and thus metastability. Indeed, 
in the model described in Ref. [12] it can be shown that 
the degree of frozen-in disorder is determined by the ra¬ 
tio of EaI^E (at fixed z/q), with increasing occupancy of 
the minor conformation below Tg with increasing ratio. 
In this regard it is interesting to note that only those 
salts found to have EaI^E above an empirical thresh- 
hold (here found computationally to be ^ 5) display an 
observable glass transition. 


V. SUMMARY AND CONCLUSION 


in Fig. [7];b), those in the A-sublattice possess close Z 

• • • H interactions with three neighbouring terminal lig¬ 
ands, while the B-sublattice EEGs are situated closer to 
the more rigid backbone Y ligands, and have close con¬ 
tacts with only one terminal SCN or N(CN )2 group. In 
the case of the non-disordered /^-(CN) 2 N(CN) 2 , the A- 
and B-sublattice display large preference for the E and 
S conformation, respectively, with 2AE/kB = 4530 and 
2010 K. For the former A-sublattice, strong coupling be¬ 
tween the EEG and thermally mobile Z = N(CN )2 lig¬ 
ands results in significant thermal displacement param¬ 
eters for the EEG, but no disorder arises, likely due to 
the strong stabilization of the eclipsed conformation. For 
the same sublattice in /i:-CuNCS, strong van der Waals 
coupling to nearby Z = SCN ligands reverses this pref¬ 
erence, favouring instead the staggered S-conformation, 
although with a much smaller 2AE/kB = 1090 K. In this 
material, the EEGs in the B-sublattice were found to 
prefer the S conformation as well, although with a larger 
2AE/kB = 1360 K. The activation energies were also ob¬ 
served to slightly differ, being 3930 and 3790 K for the 
A- and B-sublattices, respectively. Thus, the complex 
two-stage glass transition in /^-CuNCS may be related to 
subsequent freezing of the two EEG sublattices existing 
in distinct local environments. In this scenario, close Z 

• • • H contacts may play a dual role, both mediating in¬ 
teractions within the A-sublattice, and linking the two 
sublattices together as shown in Fig. [^b). A strong cou¬ 
pling of the EEG and terminal ligand vibrational degrees 


In summary, we present a comprehensive study of the 
glassy structural ordering of the ET molecules’ ethylene 
endgroups in /^-(ET) 2 X with polymeric anions X, a phe¬ 
nomenon which has important consequences for the prop¬ 
erties of the strongly-correlated electron system. Due to 
the coupling of the thermal EEG vibrations to resistance 
fluctuations, the low-frequency dynamics of the glassy 
relaxation has been determined yielding the activation 
energy of the process in samples of k-Bt and k-C\ in 
a model of non-exponential kinetics. The characteristic 
Vogel-Fulcher-Tammann behavior of the relaxation re¬ 
veals that the present materials are fragile rather than 
strong orientational glass formers. Ab initio calculations 
emphasize that the coupling between the anion layer and 
the EEG degrees of freedom plays the decisive role if a 
glass transition occurs and that the collective motion of 
EEG and anions freezes out at Tg . The estimated param¬ 
eters for the activation energies and energy differences 
between the E and S configurations explain why a dou¬ 
ble transition occurs for k-CuNCS but no transition for 
k-Cu 2 {CN)s. Predictions for the occurrence of glass-like 
transitions in certain (ET) 2 X materials are now feasible. 
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